\(\int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx\) [240]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [B] (warning: unable to verify)
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [F]
   Giac [F(-2)]
   Mupad [F(-1)]

Optimal result

Integrand size = 33, antiderivative size = 122 \[ \int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\frac {\sqrt {2} \arctan \left (\frac {\sqrt {a} \tan (e+f x)}{\sqrt {2} \sqrt {a+a \sec (e+f x)}}\right )}{\sqrt {a} (c-d) f}-\frac {2 \sqrt {d} \arctan \left (\frac {\sqrt {a} \sqrt {d} \tan (e+f x)}{\sqrt {c+d} \sqrt {a+a \sec (e+f x)}}\right )}{\sqrt {a} (c-d) \sqrt {c+d} f} \]

[Out]

arctan(1/2*a^(1/2)*tan(f*x+e)*2^(1/2)/(a+a*sec(f*x+e))^(1/2))*2^(1/2)/(c-d)/f/a^(1/2)-2*arctan(a^(1/2)*d^(1/2)
*tan(f*x+e)/(c+d)^(1/2)/(a+a*sec(f*x+e))^(1/2))*d^(1/2)/(c-d)/f/a^(1/2)/(c+d)^(1/2)

Rubi [A] (verified)

Time = 0.36 (sec) , antiderivative size = 122, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.152, Rules used = {4057, 3880, 209, 4052, 211} \[ \int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\frac {\sqrt {2} \arctan \left (\frac {\sqrt {a} \tan (e+f x)}{\sqrt {2} \sqrt {a \sec (e+f x)+a}}\right )}{\sqrt {a} f (c-d)}-\frac {2 \sqrt {d} \arctan \left (\frac {\sqrt {a} \sqrt {d} \tan (e+f x)}{\sqrt {c+d} \sqrt {a \sec (e+f x)+a}}\right )}{\sqrt {a} f (c-d) \sqrt {c+d}} \]

[In]

Int[Sec[e + f*x]/(Sqrt[a + a*Sec[e + f*x]]*(c + d*Sec[e + f*x])),x]

[Out]

(Sqrt[2]*ArcTan[(Sqrt[a]*Tan[e + f*x])/(Sqrt[2]*Sqrt[a + a*Sec[e + f*x]])])/(Sqrt[a]*(c - d)*f) - (2*Sqrt[d]*A
rcTan[(Sqrt[a]*Sqrt[d]*Tan[e + f*x])/(Sqrt[c + d]*Sqrt[a + a*Sec[e + f*x]])])/(Sqrt[a]*(c - d)*Sqrt[c + d]*f)

Rule 209

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[b, 2]))*ArcTan[Rt[b, 2]*(x/Rt[a, 2])], x] /;
 FreeQ[{a, b}, x] && PosQ[a/b] && (GtQ[a, 0] || GtQ[b, 0])

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 3880

Int[csc[(e_.) + (f_.)*(x_)]/Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)], x_Symbol] :> Dist[-2/f, Subst[Int[1/(2
*a + x^2), x], x, b*(Cot[e + f*x]/Sqrt[a + b*Csc[e + f*x]])], x] /; FreeQ[{a, b, e, f}, x] && EqQ[a^2 - b^2, 0
]

Rule 4052

Int[(csc[(e_.) + (f_.)*(x_)]*Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)])/(csc[(e_.) + (f_.)*(x_)]*(d_.) + (c_)
), x_Symbol] :> Dist[-2*(b/f), Subst[Int[1/(b*c + a*d + d*x^2), x], x, b*(Cot[e + f*x]/Sqrt[a + b*Csc[e + f*x]
])], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0]

Rule 4057

Int[csc[(e_.) + (f_.)*(x_)]/(Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)]*(csc[(e_.) + (f_.)*(x_)]*(d_.) + (c_))
), x_Symbol] :> Dist[b/(b*c - a*d), Int[Csc[e + f*x]/Sqrt[a + b*Csc[e + f*x]], x], x] - Dist[d/(b*c - a*d), In
t[Csc[e + f*x]*(Sqrt[a + b*Csc[e + f*x]]/(c + d*Csc[e + f*x])), x], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[
b*c - a*d, 0] && (EqQ[a^2 - b^2, 0] || EqQ[c^2 - d^2, 0])

Rubi steps \begin{align*} \text {integral}& = \frac {\int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)}} \, dx}{c-d}-\frac {d \int \frac {\sec (e+f x) \sqrt {a+a \sec (e+f x)}}{c+d \sec (e+f x)} \, dx}{a (c-d)} \\ & = -\frac {2 \text {Subst}\left (\int \frac {1}{2 a+x^2} \, dx,x,-\frac {a \tan (e+f x)}{\sqrt {a+a \sec (e+f x)}}\right )}{(c-d) f}+\frac {(2 d) \text {Subst}\left (\int \frac {1}{a c+a d+d x^2} \, dx,x,-\frac {a \tan (e+f x)}{\sqrt {a+a \sec (e+f x)}}\right )}{(c-d) f} \\ & = \frac {\sqrt {2} \arctan \left (\frac {\sqrt {a} \tan (e+f x)}{\sqrt {2} \sqrt {a+a \sec (e+f x)}}\right )}{\sqrt {a} (c-d) f}-\frac {2 \sqrt {d} \arctan \left (\frac {\sqrt {a} \sqrt {d} \tan (e+f x)}{\sqrt {c+d} \sqrt {a+a \sec (e+f x)}}\right )}{\sqrt {a} (c-d) \sqrt {c+d} f} \\ \end{align*}

Mathematica [A] (verified)

Time = 1.43 (sec) , antiderivative size = 144, normalized size of antiderivative = 1.18 \[ \int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\frac {2 \left (\sqrt {-c-d} \arcsin \left (\tan \left (\frac {1}{2} (e+f x)\right )\right )+\sqrt {2} \sqrt {d} \text {arctanh}\left (\frac {\sqrt {d} \tan \left (\frac {1}{2} (e+f x)\right )}{\sqrt {-c-d} \sqrt {\frac {\cos (e+f x)}{1+\cos (e+f x)}}}\right )\right )}{\sqrt {-c-d} (c-d) f \sqrt {\cos (e+f x) \sec ^2\left (\frac {1}{2} (e+f x)\right )} \sqrt {a (1+\sec (e+f x))}} \]

[In]

Integrate[Sec[e + f*x]/(Sqrt[a + a*Sec[e + f*x]]*(c + d*Sec[e + f*x])),x]

[Out]

(2*(Sqrt[-c - d]*ArcSin[Tan[(e + f*x)/2]] + Sqrt[2]*Sqrt[d]*ArcTanh[(Sqrt[d]*Tan[(e + f*x)/2])/(Sqrt[-c - d]*S
qrt[Cos[e + f*x]/(1 + Cos[e + f*x])])]))/(Sqrt[-c - d]*(c - d)*f*Sqrt[Cos[e + f*x]*Sec[(e + f*x)/2]^2]*Sqrt[a*
(1 + Sec[e + f*x])])

Maple [B] (warning: unable to verify)

Leaf count of result is larger than twice the leaf count of optimal. \(503\) vs. \(2(99)=198\).

Time = 19.78 (sec) , antiderivative size = 504, normalized size of antiderivative = 4.13

method result size
default \(\frac {\left (2 \sqrt {\left (c +d \right ) \left (c -d \right )}\, \ln \left (\csc \left (f x +e \right )-\cot \left (f x +e \right )+\sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\right ) \sqrt {\frac {d}{c -d}}+d \sqrt {2}\, \ln \left (-\frac {2 \left (-\sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, \sqrt {2}\, \sqrt {\frac {d}{c -d}}\, c +\sqrt {2}\, \sqrt {\frac {d}{c -d}}\, \sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, d +\sqrt {\left (c +d \right ) \left (c -d \right )}\, \left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right )+c -d \right )}{c \left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right )-\left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right ) d +\sqrt {\left (c +d \right ) \left (c -d \right )}}\right )-d \sqrt {2}\, \ln \left (-\frac {2 \left (\sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, \sqrt {2}\, \sqrt {\frac {d}{c -d}}\, c -\sqrt {2}\, \sqrt {\frac {d}{c -d}}\, \sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, d +\sqrt {\left (c +d \right ) \left (c -d \right )}\, \left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right )-c +d \right )}{-c \left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right )+\left (-\cot \left (f x +e \right )+\csc \left (f x +e \right )\right ) d +\sqrt {\left (c +d \right ) \left (c -d \right )}}\right )\right ) \sqrt {\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}\, \sqrt {-\frac {2 a}{\left (1-\cos \left (f x +e \right )\right )^{2} \csc \left (f x +e \right )^{2}-1}}}{2 f \sqrt {\frac {d}{c -d}}\, \left (c -d \right ) \sqrt {\left (c +d \right ) \left (c -d \right )}\, a}\) \(504\)

[In]

int(sec(f*x+e)/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))^(1/2),x,method=_RETURNVERBOSE)

[Out]

1/2/f/(d/(c-d))^(1/2)/(c-d)/((c+d)*(c-d))^(1/2)/a*(2*((c+d)*(c-d))^(1/2)*ln(csc(f*x+e)-cot(f*x+e)+((1-cos(f*x+
e))^2*csc(f*x+e)^2-1)^(1/2))*(d/(c-d))^(1/2)+d*2^(1/2)*ln(-2*(-((1-cos(f*x+e))^2*csc(f*x+e)^2-1)^(1/2)*2^(1/2)
*(d/(c-d))^(1/2)*c+2^(1/2)*(d/(c-d))^(1/2)*((1-cos(f*x+e))^2*csc(f*x+e)^2-1)^(1/2)*d+((c+d)*(c-d))^(1/2)*(-cot
(f*x+e)+csc(f*x+e))+c-d)/(c*(-cot(f*x+e)+csc(f*x+e))-(-cot(f*x+e)+csc(f*x+e))*d+((c+d)*(c-d))^(1/2)))-d*2^(1/2
)*ln(-2*(((1-cos(f*x+e))^2*csc(f*x+e)^2-1)^(1/2)*2^(1/2)*(d/(c-d))^(1/2)*c-2^(1/2)*(d/(c-d))^(1/2)*((1-cos(f*x
+e))^2*csc(f*x+e)^2-1)^(1/2)*d+((c+d)*(c-d))^(1/2)*(-cot(f*x+e)+csc(f*x+e))-c+d)/(-c*(-cot(f*x+e)+csc(f*x+e))+
(-cot(f*x+e)+csc(f*x+e))*d+((c+d)*(c-d))^(1/2))))*((1-cos(f*x+e))^2*csc(f*x+e)^2-1)^(1/2)*(-2*a/((1-cos(f*x+e)
)^2*csc(f*x+e)^2-1))^(1/2)

Fricas [A] (verification not implemented)

none

Time = 0.49 (sec) , antiderivative size = 963, normalized size of antiderivative = 7.89 \[ \int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\text {Too large to display} \]

[In]

integrate(sec(f*x+e)/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))^(1/2),x, algorithm="fricas")

[Out]

[-1/2*(sqrt(2)*sqrt(-1/a)*log((2*sqrt(2)*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*sqrt(-1/a)*cos(f*x + e)*sin(f
*x + e) + 3*cos(f*x + e)^2 + 2*cos(f*x + e) - 1)/(cos(f*x + e)^2 + 2*cos(f*x + e) + 1)) + sqrt(-d/(a*c + a*d))
*log(-((c^2 + 8*c*d + 8*d^2)*cos(f*x + e)^3 + (c^2 + 2*c*d)*cos(f*x + e)^2 - 4*((c^2 + 3*c*d + 2*d^2)*cos(f*x
+ e)^2 - (c*d + d^2)*cos(f*x + e))*sqrt(-d/(a*c + a*d))*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*sin(f*x + e) +
 d^2 - (6*c*d + 7*d^2)*cos(f*x + e))/(c^2*cos(f*x + e)^3 + (c^2 + 2*c*d)*cos(f*x + e)^2 + d^2 + (2*c*d + d^2)*
cos(f*x + e))))/((c - d)*f), -1/2*(sqrt(2)*sqrt(-1/a)*log((2*sqrt(2)*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*s
qrt(-1/a)*cos(f*x + e)*sin(f*x + e) + 3*cos(f*x + e)^2 + 2*cos(f*x + e) - 1)/(cos(f*x + e)^2 + 2*cos(f*x + e)
+ 1)) + 2*sqrt(d/(a*c + a*d))*arctan(2*(c + d)*sqrt(d/(a*c + a*d))*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*cos
(f*x + e)*sin(f*x + e)/((c + 2*d)*cos(f*x + e)^2 + (c + d)*cos(f*x + e) - d)))/((c - d)*f), -1/2*(sqrt(-d/(a*c
 + a*d))*log(-((c^2 + 8*c*d + 8*d^2)*cos(f*x + e)^3 + (c^2 + 2*c*d)*cos(f*x + e)^2 - 4*((c^2 + 3*c*d + 2*d^2)*
cos(f*x + e)^2 - (c*d + d^2)*cos(f*x + e))*sqrt(-d/(a*c + a*d))*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*sin(f*
x + e) + d^2 - (6*c*d + 7*d^2)*cos(f*x + e))/(c^2*cos(f*x + e)^3 + (c^2 + 2*c*d)*cos(f*x + e)^2 + d^2 + (2*c*d
 + d^2)*cos(f*x + e))) + 2*sqrt(2)*arctan(sqrt(2)*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*cos(f*x + e)/(sqrt(a
)*sin(f*x + e)))/sqrt(a))/((c - d)*f), -(sqrt(d/(a*c + a*d))*arctan(2*(c + d)*sqrt(d/(a*c + a*d))*sqrt((a*cos(
f*x + e) + a)/cos(f*x + e))*cos(f*x + e)*sin(f*x + e)/((c + 2*d)*cos(f*x + e)^2 + (c + d)*cos(f*x + e) - d)) +
 sqrt(2)*arctan(sqrt(2)*sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*cos(f*x + e)/(sqrt(a)*sin(f*x + e)))/sqrt(a))/
((c - d)*f)]

Sympy [F]

\[ \int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\int \frac {\sec {\left (e + f x \right )}}{\sqrt {a \left (\sec {\left (e + f x \right )} + 1\right )} \left (c + d \sec {\left (e + f x \right )}\right )}\, dx \]

[In]

integrate(sec(f*x+e)/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))**(1/2),x)

[Out]

Integral(sec(e + f*x)/(sqrt(a*(sec(e + f*x) + 1))*(c + d*sec(e + f*x))), x)

Maxima [F]

\[ \int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\int { \frac {\sec \left (f x + e\right )}{\sqrt {a \sec \left (f x + e\right ) + a} {\left (d \sec \left (f x + e\right ) + c\right )}} \,d x } \]

[In]

integrate(sec(f*x+e)/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))^(1/2),x, algorithm="maxima")

[Out]

integrate(sec(f*x + e)/(sqrt(a*sec(f*x + e) + a)*(d*sec(f*x + e) + c)), x)

Giac [F(-2)]

Exception generated. \[ \int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\text {Exception raised: TypeError} \]

[In]

integrate(sec(f*x+e)/(c+d*sec(f*x+e))/(a+a*sec(f*x+e))^(1/2),x, algorithm="giac")

[Out]

Exception raised: TypeError >> an error occurred running a Giac command:INPUT:sage2:=int(sage0,sageVARx):;OUTP
UT:index.cc index_m i_lex_is_greater Error: Bad Argument Value

Mupad [F(-1)]

Timed out. \[ \int \frac {\sec (e+f x)}{\sqrt {a+a \sec (e+f x)} (c+d \sec (e+f x))} \, dx=\int \frac {1}{\cos \left (e+f\,x\right )\,\sqrt {a+\frac {a}{\cos \left (e+f\,x\right )}}\,\left (c+\frac {d}{\cos \left (e+f\,x\right )}\right )} \,d x \]

[In]

int(1/(cos(e + f*x)*(a + a/cos(e + f*x))^(1/2)*(c + d/cos(e + f*x))),x)

[Out]

int(1/(cos(e + f*x)*(a + a/cos(e + f*x))^(1/2)*(c + d/cos(e + f*x))), x)